arXiv:1507.07550v2 [hep-ph] 14 Dec 2015 


FTU AM-15-15 IFT-UAM/CSIC-15-057 LPT-Orsay-15-58 
Prepared for submission to JHEP 


On the Viability of Minimal Neutrinophilic 
Two-Higgs-Doublet Models 


P. A. N. Machado,® Y. F. Perez , 6 O. Sumensari, 6,c Z. Tabrizi, d,e and R. Zukanovich 
Funchal 6 

“Departamento de Fisica Tedrica and Instituto de Fisica Teorica, IFT-UAM/CSIC, 

Universidad Autonoma de Madrid, Cantoblanco, 28049, Madrid, Spain 
b Departamento de Fisica Matemdtica, Instituto de Fisica, Universidade de Sao Paulo, 

C. P. 66.318, 05315-910 Sao Paulo, Brazil 
c Laboratoire de Physique Theorique (Bat.210), 

Universite Paris Sud and CNRS (UMR 8621), F-91405 Orsay-Cedex, France 
d School of Particles and Accelerators, 

Institute for Research in Fundamental Sciences (IPM), P.O.Box 19395-1195, Tehran, Iran 
e Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas (UNICAMP), Rua Sergio 
Buarque de Holanda, 111, Campinas, SP, 13083-859, Brazil 

E-mail: pedro.machado@uam.es, yfperezg@if.usp.br, 

olcyr.sumensari@usp.br, tabrizi,physics@ipm.ir, zukanov@if.usp.br 

Abstract: We study the constraints that electroweak precision data can impose, after the 
discovery of the Higgs boson by the LHC, on neutrinophilic two-Higgs-doublet models which 
comprise one extra SU( 2) x 1/(1) doublet and a new symmetry, namely a spontaneously 
broken Z 2 or a softly broken global t/(l). In these models the extra Higgs doublet, via 
its very small vacuum expectation value, is the sole responsible for neutrino masses. We 
find that the model with a Z 2 symmetry is basically ruled out by electroweak precision 
data, even if the model is slightly extended to include extra right-handed neutrinos, due 
to the presence of a very light scalar. While the other model is still perfectly viable, the 
parameter space is considerably constrained by current data, specially by the T parameter. 
In particular, the new charged and neutral scalars must have very similar masses. 
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1 Introduction 

The smallness of neutrino masses suggests a mass generating mechanism distinct from the 
usual Higgs mechanism, which resides in a scale different from the electroweak one. From 
neutrino oscillation experiments, we know that neutrinos are massive and that mass and 
flavor eigenstates do not coincide. Besides, other terrestrial [1, 2] and cosmological [3, 4] 
experiments indicate that neutrino masses should be below the eV scale. Therefore, if the 
same Higgs mechanism is responsible for the top and neutrino masses, then the Yukawa 
couplings would span twelve orders of magnitude, evincing an unpleasant and inexplicable 
hierarchy. 

A well known alternative is the seesaw mechanism [5-7]. In this scenario, the light 
neutrino masses are suppressed by some heavy physics, for instance, right-handed Majorana 
neutrino masses [5, 6, 8, 9]. What typically happens is that the scale at which new physics 
can be found is extremely high, much above the TeV scale, rendering the model intangible, 
except for the possible presence of neutrinoless double beta decay 1 . The latter could 
also originate from some physics that do not comprise the main contribution to neutrino 
masses [13, 14], and hence it does not consist of a test of the seesaw mechanism by itself. 

Another possibility is to generate neutrino masses by a copy of the Higgs mechanism, 
having a second Higgs doublet, but with a much smaller vacuum expectation value (vev). 
This can be achieved in a two-Higgs-doublet model (2HDM) where one of the scalars gives 

Nevertheless, there are alternative models which exhibit a low scale, as for instance the inverse seesaw 
scenario [10-12]. 
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mass to the charged fermions, while the other one acquires a very small vev and generates 
neutrino masses with 0(1) Yukawa couplings, a neutrinophilic 2HDM. As a consequence, 
neutrino masses would generically require new physics at the TeV scale (or even lower). For 
instance, by imposing a lepton number symmetry and adding three right-handed neutrinos 
which carry no lepton number, a type I seesaw mechanism can be realized below the TeV 
scale [15]. Moreover, lepton number could be conserved and a Z 2 symmetry [16, 17] or 
a global U( 1) [18] could be used to prevent the SM Higgs boson to couple to neutrinos, 
yielding Dirac neutrinos. Also, the 2HDM could be augmented by a type III seesaw and a 
fj, — t symmetry, giving rise to interesting LHC phenomenology [19]; or by a singlet scalar 
and a Z 3 symmetry, possibly generating lepton flavor violating signals [20] . It is important 
to note that such models are stable against radiative corrections [ 21 , 22 ]. 

On general grounds, a new symmetry is typically invoked to prevent the first scalar 
doublet from coupling to neutrinos as well as to enforce the second one to interact only 
with them. These models introduce a minimal new field content which should materialize 
as particles below the TeV scale. The presence of such a low scale in the theory might 
have important phenomenological consequences, like the presence of light scalar particles 
(for instance, supernova energy loss strongly constrains such scenarios [23]). After the 
discovery of a 125 GeV scalar by the LHC experiments, new limits from electroweak pre¬ 
cision data can be derived on the allowed parameter space of such models. The purpose of 
this manuscript is to investigate to what extent these minimal neutrinophilic 2HDMs can 
survive electroweak precision data scrutiny. 

In sec. 2 we briefly review the neutrinophilic 2HDMs which we will study in this work. 
In sec. 3 we describe the theoretical and experimental constraints that will be imposed on 
these models in sec. 4. Finally, in sec. 5 we present our conclusions. 


2 Neutrinophilic Two-Higgs-Doublet Models 


We first start by making general considerations on the 2HDM and the link to neutrino 
masses. The most general scalar potential for a 2HDM is 


V (< 3 ?i, <F 2 ) = 1 + m| 2 d>^<I>2 - (m\ 2 <&\$2 + h.c.) 

+ y (<f{<Fi) 2 + y ( 4 > 5 <F 2 ) 2 + A3<FI<L1<F5<L2 + A4<FI<F2<F2<F1 


+ 


y($I*2) 2 + 


<f|<Fi + A 7 4>^4>2 ) <Fj<F 2 + h.c. 


st 


st , 


( 2 . 1 ) 


where <Fi and <F 2 are two scalar doublets with hypercharge Y = +1. For the vacuum 
expectation values of the two scalars, we adopt the notation (<Fi) = v\/y/2, (<F 2 ) = V2/V2, 
and we pick <F 2 to be the one responsible for neutrino masses. In order to have sizable 
Yukawa coupling for neutrinos, it is required that u 2 -C v\ ~ 246 GeV = v, where v 2 = 
v 2 + v 2 ■ In principle, the parameters m 2 2 , A 5 , Ag, and A 7 can be complex. Nevertheless, in 
all models we analyze, the symmetries will forbid both Ag and A 7 , and only m\ 2 or A 5 will 
be allowed to be non-zero. A single phase of the aforementioned parameters can always be 
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absorbed in a redefinition of the scalar fields, and therefore we can take all scalar potential 
parameters to be real without loss of generality. 

To forbid the coupling between neutrinos and <f>i, a symmetry is called for. In 
this minimal setup, there are two straightforward examples. The first possibility is a 
Z 2 symmetry under which only <I >2 and the right-handed neutrinos are charged, forcing 
mi 2 = A 6 = A 7 = 0. An alternative is to trade the Z 2 by a global 17(1), yielding, in princi¬ 
ple, mi 2 = A 5 = A 6 = A 7 = 0. In this case, to avoid the presence of a massless Goldstone 
boson, a soft breaking is introduced by having a non-zero but small mi 2 - On the other 
hand, if a softly broken Zjv, N > 2, symmetry is postulated, A 5 might be forcefully zero 
as well, making this case identical to the 17(1) scenario. Therefore, the phenomenology 
of a softly broken Zj\r >2 model is identical to the softly broken U (1) case. Anyhow, in 
all realizations we will study here Ag = A 7 = 0 , so these couplings will be disregarded 
henceforth. 

One last option that one could consider would be to gauge the U (1) symmetry, avoiding 
the massless Goldstone boson. Nevertheless, in such a scenario, the corresponding gauge 
boson as well as one of the neutral scalars would be extremely light, with mass around the 
V 2 scale. This seems, at first glance, phenomenologically quite problematic. We do not 
investigate this possibility here as it would require a completely different study compared 
to the other two cases. 

The two complex scalar SU (2) doublets can be written as 


$a = 


(v a + Pa + iria)/V 2 , 


a = 1 , 2 . 


( 2 . 2 ) 


After electroweak symmetry breaking, three Goldstone bosons become the longitudinal 
modes of the W and Z bosons. Then, the remaining scalar spectrum is composed of two 
charged particles, i7 ± , two CP-even neutral bosons, h and H, and one CP-odd neutral 
boson, A. The physical fields are given by 

H + = (f)^ sin (3 — <f >2 cos (3, A = rji sin (3 — rj 2 cos /3, (2.3) 

h = — pi cos a — P 2 sin a, H = p\ sin a — P 2 cos a, (2.4) 


where the angles a and (3 are associated with the rotations that diagonalize the mass 
matrices 


tan( 2 a) 
tan f3 


_ 2(~mf 2 + A 345 V 1 V 2 ) _ 

m 2 l 2 (v 2 /vi - V1/V2) + Ait>j - \ 2 vl ’ 

V2 

Vi 


(2.5) 

( 2 . 6 ) 


where A 345 = A 3 + A 4 + A 5 . We will see below that both a and f3 are expected to be very 
small. Hence, h behaves very similarly to the SM Higgs, while the scalars H, A, develop 
neutrinophilic interactions in the Yukawa sector, as described below 


£ y = 1 —Hu i u i - i— Ava^vi - ^ mvi [U* ti H + UiP L l + h.c.], 

V 2 V 2 V 2 


(2.7) 
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where m Vi are neutrino masses and Ua is the PMNS matrix. As we will see in sec. 4, 
the tree-level stationary conditions on the potential, dV/d&i = 0 , can be used to write 
the diagonal mass parameters mu as functions of m 2 2 , the quartic couplings and the vevs. 
With that in mind, we can consider the quartic couplings as free parameters and express 
them in terms of the physical masses, vevs and mixing angles [24] (see appendix A). 

Next we describe the two specific realizations of the neutrinophilic scenarios that will 
be studied in this paper. 

2.1 Neutrinophilic 2HDM: Z 2 symmetry 

The model to be studied was proposed by Gabriel and Nandi [16] 2 . It consists of a 2HDM 
where both the right-handed neutrinos and one of the scalar doublets, 4 > 2 , are charged under 
a Z 2 symmetry. The consequence is that the masses of the charged fermions come solely 
from the $1 vev, and neutrinos, which are Dirac fermions in this scenario as the authors 
impose lepton number conservation, couple exclusively to <f> 2 . This extra symmetry can, in 
principle, be dropped allowing for Majorana neutrinos with a low scale realization of the 
seesaw mechanism. We will also investigate this possibility in our analysis. 

In the scalar potential (2.1) of this model, the parameters m\ 2 and A 6,7 will vanish 
due to the Z 2 symmetry. The smallness of neutrino masses is explained by the very low 
scale at which Z 2 is broken, preferably V 2 < 0(eV ). 3 A tiny vifv\ ratio and the absence 
of an explicit breaking m 2 2 term leads to almost no mixing between the doublets. The 
smallness of tan/3 and tana can be seen from eqs. (2.5) and (2.6) after imposing V 2 /V 1 —> 0. 
Therefore, apart from its couplings to neutrinos, 4>i behaves almost identically to the SM 
Higgs doublet, so we do not expect any observable deviation from the Higgs couplings to 
the SM particles, except possibly the loop induced couplings, e.g. hy y. 

The second doublet displays some interesting features. Through the Yukawa coupling, 
the neutral components couple almost only to neutrinos, while the charged scalars mediate 
interactions between neutrinos and charged leptons (see eq. (2.7)). The Yukawas are ideally 
expected to be of 0(1). The neutral scalars couple to the W and Z bosons, but notice 
that triple gauge couplings (TGCs) involving only one scalar are highly suppressed by the 
small vev, V 2 - Obviously, TGCs with two scalars and one gauge boson are present and may 
provide a sizeable pair production cross section at colliders, for instance pp —> A* —> H + H~ 
at the LHC. 

The scalar spectrum of this model is quite constrained. By setting m\ 2 = 0 in 
eqs. (A. 6 )-(A. 8 ), as well as sin 2 a, sin 2 (3 <C 1, we notice that: (i) h is identified as the 
125 GeV Higgs particle found at the LHC, and this essentially fixes Ai ~ 0.26 (see eq. (A.l)); 
(ii) the neutrinophilic neutral scalar H is extremely light, mu ~ 0(v 2 ) -C v\ and (iii) for 
not so large values of the quartic couplings, the charged scalars and the pseudoscalar masses 

2 The same model was previously also discussed in ref. [25] where the focus was on the origin of the 
second doublet from neutrino condensation. 

3 It is known that the breaking of discrete symmetries leads to the formation of domain walls, which may 
store unacceptably large quantities of energy, unless the vev responsible for this breaking is below d?(10 -2 ) 
GeV [26, 27]. Nevertheless, as the second scalar has to have a vev small enough to explain neutrino masses, 
domain walls do not pose a bound on neutrinophilic 2HDMs. 
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are bounded to be about or below the TeV scale. When we analyse the viability of this 
model in sec. 4, it will turn out that oblique parameters will play a decisive role in con¬ 
straining it, due to the peculiar structure of the scalar spectrum. The sensitivity of the S 
parameter to the presence of a very light neutral scalar, m# ~ 0(v 2 ), will essentially rule 
out the model. 


2.2 Neutrinophilic 2HDM: softly broken global 17(1) symmetry 

The second model we study was proposed by Davidson and Logan [18]. Analogously to 
the other scenario, both $2 and right-handed neutrinos are charged under a new global 
17(1). The model spans = 0 and a small m 2 2 which breaks the symmetry softly and 
generates neutrino masses. The presence of the soft breaking mass term, is required in 
order to avoid a massless Goldstone boson which might create problems with cosmology 
and electroweak precision data. Neutrinos are Dirac particles, as the Majorana mass term 
is strictly forbidden by the new 17(1). From eq. (A. 8), we write 

m i2 = s i n P cos P (2-8) 


and we observe that to obtain simultaneously V 2 ~ eV and tha ~ 0(100 GeV) one would 
need m\ 2 ~ (200 keV) 2 . As said before, to avoid the issues of having a massless Goldstone, 
instead of softly breaking the new 17(1) symmetry, one could also envisage to gauge it. 
Nonetheless, the theory would contain a very light vector resonance as a consequence of the 
small vev, and it is not clear if such a model can satisfy all neutrino data and astrophysical 
constraints. We do not explore this possibility here. 

The presence of a non-zero m\ 2 term makes this case fairly different from the last one. 
From eq. (A.7), we notice that the mass of the neutrinophilic scalar, mn, increases with 
M, and therefore the H mass in this scenario is not bounded by V 2 as in the previous case. 
As we will see later, this will ease the constraints from the oblique parameters. Combining 
eq. (A. 2) with the definition M 2 = m 2 2 /(sin/3 cos/3), and imposing tan/3 = V 2 /V -C 1, we 
obtain 


A2 = 


cot 2 /3 M 2 + 


cos 2 a 
sin 2 /3 


m H + 


sin 2 a 
sin 2 /3 



1 f 2 1 v \ 

- ~2 rn H -m 12 — 
V 2 \ V 2 ) 


+ 


sin 2 a 
sin 2 /3 v 2 ’ 


(2.9) 


which indicates that 

\ m H ~ ml 2 v/v 2 \ < 0(v 2 ). (2.10) 

To grasp the impact of this conclusion, assume that ?b 2 2 = m 2 H v 2 /v. Hence, from eq. (A.8) 
we see that nriA ~ run, so the neutrinophilic CP-odd and CP-even scalars, A and H, are 
degenerate in mass. We emphasize that this degeneracy by itself is not a fine tuning of 
the model: the degenerate spectrum arises naturally given the symmetries of the scalar 
potential and the hierarchy between the vevs. As a last comment, we emphasize that since 
m 2 2 is the only source of f7(l) breaking, it is natural in the t’Hooft sense - m 2 2 only receives 
radiative corrections proportional to itself [21, 22]. 
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3 Theoretical and Experimental Electroweak Data Constraints 
3.1 Theoretical Constraints 

There are a number of conditions to be fulfilled by the scalar potential. These will be used 
to constrain the parameter space, ultimately restricting the range of physical scalar masses, 
having an important impact on the phenomenology of the models. To have stability at 
tree-level, the following constraints should be fulfilled [28] 

Ai ,2 > 0, A 3 > — (A 1 A 2) 1 / 2 and A 3 + A 4 — |A 5 1 > — (A 1 A 2 ) 1 / 2 . (3.1) 


In addition, the stationary conditions dV/d^i = 0 read 

+ mhv 1 - m\ 2 v 2 = 0 , 

2 2 (3.2) 

^2 3 ^345 2,2 2 n 

—n 2 + -^-v 2 v 1 + m 22 v 2 - m 12 vi = 0, 

which allow us to write m- as functions of m \ 2 , v\ and v 2 . If m\ 2 = 0, it is easy to see 
that there are at least two equivalent stable solutions, (v , 0 ) or ( 0 , v ) (although they may 
not be the global minima). In this case, the vev is precisely the electroweak scale, one of 
the scalars is exactly the Higgs and the other one is inert. For m\ 2 7 / 0, these equations 
cannot be solved analytically. Nevertheless, if rn\ 2 <C v 2 a perturbative approach yields 


v\ ~ v, 


m■■ 


v 2 


12 


-V345 „,2 


v z + m 


-v, 


22 


(3.3) 


and a symmetric solution interchanging the indices 1 -H- 2 , which reveals that the small vev 
necessary to satisfactorily explain small neutrino masses might require a correspondingly 
small m\ 2 parameter. This can be understood intuitively, as the breaking of the 1/(1) 
happens only through the soft breaking term m\ 2 . In general, there can be more than one 
solution satisfying the stationary conditions (3.2), and hence different non-trivial and non¬ 
degenerate minima (vi,v 2 ) and (v \, v' 2 ) might coexist. It is possible to check analytically 
if the chosen vacuum is the deepest one in the potential for a 2HDM with = 0 [29]. 
In this case, the potential describes a Z 2 symmetry softly broken by m\ 2 . Both models we 
deal with here are special cases of such scenario. In the absence of an explicit breaking, 
that is m\ 2 = 0 , there can be multiple minima, but they are degenerate and hence stability 
is not threatened. This is the case of the h 2 model we analyze. For the softly broken 1/(1) 
model, it can be shown that the chosen vacuum is the deepest one (at tree-level) if and 
only if the following condition is satisfied [29]: 


D = mf 2 (?nf 1 — K 2 m, 22 )( tan (5 — k) > 0, (3.4) 

with n = ^/\\/\ 2 . Although for a general 2HDM scenario this bound may be important, 
for the neutrinophilic case we have checked that it does not lead to any significant effect 
on the parameter space, after the other constraints are taken into account, but we include 
it in the analysis of the softly broken U{1) model for completeness. 
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Another theoretical requirement is to satisfy the tree-level pertubative unitarity con¬ 
dition [30-32]. If the quartic couplings are too large, the lowest order amplitudes for 
scalar-scalar scattering may violate unitarity at high enough scales, requiring additional 
physics to mitigate this issue. To obtain the constraint, the scalar-scalar S matrix is 
computed and the following conditions are imposed on its eigenvalues 

|a±|, |&±|, |c±|, |/±|, |e lj2 |, |/i|, |pi| < 8vr, ( 3 . 5 ) 

where 


a± = -(Ai + A 2 ) ± \J -(Ai - A 2 ) 2 + (2A 3 + A4) 2 , 

( 3 . 6 a) 

b± = 2(^1 + ^2) ± 2\A^i “ ^) 2 + 

( 3 . 6 b) 

C± = t,(Ai + A 2 ) ± ^y/(Ai - A 2 ) 2 + 4 A|, 

( 3 . 6 c) 

/+ = A3 + 2A4 + 3 A. 5 , 

( 3 . 6 d) 

/- = A 3 + A5, 

( 3 . 6 e) 

ei = A3 + 2A4 — 3A.5, 

( 3 . 6 f) 

= A3 — A5, 

( 3 . 6 g) 

fi = A3 + A4, 

( 3 . 6 h) 

Pi = A3 — A4. 

( 3 . 6 i) 


To have an idea of the impact of these bounds, one can conservatively assume that all | A*| 
should be smaller than 87 r (the actual bound is always more stringent than that). Some 
authors prefer to use a stronger limit of 47 r. We checked that this does not change very 
much the allowed regions. 

Evidently, even if tree-level unitarity is satisfied, loop corrections could still play an 
important role leading to violation of unitarity at some scale and thus demanding the 
presence of new physics below such energies. This could be particularly relevant when some 
of the tree-level constraints are just barely satisfied, as the size of the quartic couplings 
could enhance the loop contributions. Nevertheless, we only take into account unitarity 
constraints at tree-level, as a full one loop evaluation of the parameter space is beyond the 
scope of this manuscript. 

3.2 Electroweak Data Constraints 

Oblique Parameters. The impact of a second Higgs doublet in the so-called electroweak 
precision tests (EWPT), encoded in the Peskin-Takeuchi parameters S, T, and U [33], has 
been studied in the literature to a great extent (see for instance refs. [34-36]). These are 
radiative corrections to the gauge boson two point functions, known as oblique corrections. 
For the precise expressions of S, T, U, we point the reader to the aforementioned references. 

The S parameter encodes the running of the neutral gauge bosons two point functions 
(ZZ, Z 7 and 77 ) between zero momentum and the Z pole. Therefore, it should be specially 
sensitive to new physics at low scales, particularly below the Z mass. Thus, we expect it to 
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be important in the presence of very light neutral scalars, as is the case for the Z 2 model. 
The T parameter measures the breaking of custodial symmetry at zero momentum, that is, 
the difference between the WW and the ZZ two point functions at q 2 = 0. It usually plays 
a significant role in constraining the parameter space of particles charged under SU(2)l. 
Splitting the masses of particles in a doublet breaks custodial symmetry and affects T. 
As we will see later, in the softly broken U( 1) scenario, the T parameter will provide the 
major constraint on the mass splitting m#± — mA-, forcing the scalar spectrum of this 
model to be somewhat degenerate. Last, and this time least, the U parameter (or better, 
the combination S ± U) is somewhat similar to S but for the W bosons, being sensitive to 
light charged particles in the loops. Given the fact that light charged particles are excluded 
by LEP data [37, 38], usually U is the least important of these three precision parameters, 
having a minor impact on the model phenomenology, we have checked that this is indeed 
the case for all scenarios analyzed here. 

To evaluate the impact of the EWPT on the neutrinophilic 2HDM scenarios, we cal¬ 
culate S, T, and U using the results available in ref. [36], and we use the latest GFITTER 
values for the best fit, uncertainties and covariance matrix [39], 


A S SM = 0.05 ±0.11, 
A T sm = 0.09 ±0.13, 
A U SM = 0.01 ±0.11, 


V = 


l 1 0.90 -0.59 

0.90 1 -0.83 

y -0.59 -0.83 1 


composing the \ 2 function as 


X 2 = 




SM\ 


hJ 


(3.7) 


(3.8) 


with Xi = AS, AT, AU and the covariance matrix < 7 ^ = (JiVijCTj, in which (( 71 ,( 72 ,( 73 ) = 
(0.11,0.13,0.11). As we are interested in the goodness of fit of the model to the EWPT 
data, the 1, 2, and 3(7 regions are calculated using y 2 = 3.5,8.0,14.2, respectively. 

Higgs invisible width. When the first doublet acquires a vev, triple scalar vertices 
like hSS ( S = H, A) are induced. Therefore, light neutral scalars with 2ms < m h. could 
contribute to the Higgs invisible width h —> SS, and sequentially S —> vv. Because of the 
small tan/3 of the model, the Higgs boson couplings to the Standard Model particles is 
basically unchanged. Hence, the contribution to the Higgs total width due to the invisible 
decay will suppress all Standard Model branching fractions by the ratio r^ M /T)] ew . In this 
scenario, as the only modification to the Higgs branching fractions is the addition of an 
invisible channel, the LHC 8 TeV data bound is BR (h —> invisible) < 0.13 at 95% CL [40]. 

In our framework, the decay rate of such a process is given by [41] 


T(h-+SS) 


9hSS 

32irm.h 



(3.9) 
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with 


9hAA = 


2v 


2 2 sin(a-3/3) 


(2 rriA-ml)- 


sin 2/3 


+ (8 mf 2 ~ sin2/3(2m^ + 3 m^)) 


sin(/3 + a) 


sin 2 2/3 


9hHH = — cos(/3 — a) 
v 


2m 


12 


sin 2/3 


+ ( 2m H + rn h - - 


6 m 


h 


12 


sin 2a 


sin 2/3 J sin 2/3 


(3.10) 

(3.11) 


While the couplings between the SM Higgs and the SM fermions, ghff = mf/v, are well 
below one due to the suppression by the EW scale (except for the top, to which the Higgs 
cannot decay), the trilinear scalar couplings are typically much larger, ghss ~ m h/ v ~ 
60GeV, unless there is some sort of cancellation happening [41]. Therefore, SM Higgs 
decays to lighter scalars may have an important phenomenological impact, see e.g. ref. [42], 
specially because the total Higgs width in the Standard Model is predicted to be very small, 
around 4.07 MeV [43] . 

Higgs to diphoton. The charged scalars will contribute to the h —> 77 width, and 
thus we also analyse the impact on this observable 1 . The h diphoton width is a destructive 
interference effect mainly between W and top loops, where the latter dominate. Charged 
scalars contribute with the same sign as the W, and their contribution usually do not 
overcome the top one. Therefore, we expect h -A- 77 to be somewhat suppressed in most 
cases. The expression for the h —> 77 width at one loop can be found in many papers, see, 
for instance ref. [44]. For reference, the current ATLAS+CMS combination value of the 
Higgs to diphoton signal strength is /r 77 = Tiding [45]. 

Z invisible width. We also have to consider possible extra contributions to the Z 
invisible width coming from the decays Z -A Svv with S = A, H and ms < mz■ In the 
model with a softly broken U{ 1) symmetry, the expression for T(Z -A Suv) = T(Z -A 
Auv) + r(Z —> Hvv) can be easily calculated and reads 


T(Z^SvP) = 


384ir 3 m 4 5 z 


9s{q 2 ) + 


2 2 

nl u,tot 


"( mz-ms ) 2 


2 COS 6\v 

fs(9 2 


d q 2 


X 1/2 (9 2 , 


m z , mt 


(q 2 — m |) 2 + w,|T| 


/ o 9 9 

-1 ( m z + — q 


X 1 / 2 (<7 2 , m%,m 2 z ) 


l v A 1 / 2 (g 2 ,m|,m|) 


(3.12) 


where ms is the mass of the neutrinophilic scalars H and A, which are degenerate in mass, 
and the total width is given by 


_ ms ml tot 

87T v\ 


(3.13) 


4 The h -A decay will also be modified, but due to the smaller branching ratio and subsequent 

suppression by requiring the Z to decay leptonically, we do not expect it to provide any significant sensitivity 
in the near future. 
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We also define A(a 2 ,& 2 ,c 2 ) = (a 2 — (6 —c) 2 )(a 2 — (6 + c) 2 ) and 

fs(q 2 ) = 4m| [(m| - g 2 )(m| - m|g 2 + g 4 + m|(m| - 4g 2 )) + + m| - g 2 )] , 

(3.14) 

9s{q 2 ) = 4m|(g 2 - m|) + m|[4m|(2g 2 - r|) + g 2 (r| - 8 g 2 )] + g 2 (m| - 8 ?n|g 2 + 4g 4 ). 

(3.15) 

The ratio between m 2 tot = ^ m 2 . and u 2 arrives from the neutrino Yukawas. Clearly, if 
the Yukawas are small, both widths vanish, so we expect this bound to be more significant 
for lower v 2 and larger neutrino masses. To constrain extra contributions from new physics 
to the Z invisible width, we use LEP result T exp (Z —> invisible) = 499.0(15) MeV and 
the Standard Model prediction T SM (Z —> invisible) = 501.69(6) MeV [38], which yields 
r NP (Z — > invisible) < 1.8 MeV at 3cr (notice that there is a mild 2 a discrepancy between 
the data and the SM predicted value). In the case of the Z 2 symmetry model, one must 
take care while doing the computation, since run mz, as the expression for the width 
has an infrared divergence, which cancels out with radiative conditions. As we will see in 
sec. 4, the other constraints will exclude most of the parameter space of this model. For 
this reason we will not discuss the constraints from the Z invisible width in this scenario. 

Collider bounds on charged scalars. The charged scalars can be pair produced 
directly at colliders via s-channel off shell photon or Z exchange. Due to the neutrinophilic 
character of the second Higgs doublet and small admixture with the SM degrees of freedom, 
the charged scalars decay almost only to Iv. Therefore, we use the corresponding LEP 
bound, i.e. m H ± > 80 GeV [37, 38]. 

It is not clear how LHC data improves the situation. There has been some studies 
on the LHC sensitivity to such charged scalars, mainly focused on 14 TeV center of mass 
energy [46-50], but to the best of our knowledge, there has been no dedicated experimental 
search for charged scalars in neutrinophilic 2HDMs. As V 2 is very small, the main pro¬ 
duction modes of H± would be pair production through vector boson fusion or off-shell 
s-channel photon and Z exchange, and the tipical t —>■ H + b would be absent due to small 
tan (3. The LHC sensitivity then would come mainly from opposite sign dilepton plus miss¬ 
ing energy, which has SM W pair production as an irreducible background. Moreover, 
the branching ratios of the charged scalar depend on the neutrino masses and the mass 
ordering. If the tv branching ratio is dominant, the sensitivity is expected to be smaller. 
Therefore, to be conservative, we will scan the parameter space considering only the LEP 
bound. 

Anomalous magnetic moments and other constraints. In principle, the charged 
scalars could also contribute to charged lepton g —2 values, but the corresponding amplitude 
at one loop is suppressed by mj/m^ ± (see ref. [51] for a recent analysis on the impact of a 
second Higgs doublet on the muon g — 2). We have checked that the 1-loop contribution to 
both muon and electron g — 2 is negligible due to that suppression, while the tau g — 2 is not 
measured with enough precision to pose a bound. For a general 2HDM, it has been noticed 
that two loop Barr-Zee diagrams [52] can be more important than 1-loop contributions, 
but this is not the case in the neutrinophilic 2HDM, as the charged lepton couplings to H 
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Figure 1. Neutrinophilic 2HDM with Z 2 symmetry. The red points are allowed by electroweak 
precision data (oblique parameters) at 3cr, while the gray points are ruled out at 3cr or more in the 
S x T plane. No point was found within the 2cr region. 

and A are suppressed by tan/3 5 * . Therefore we conclude that the electron, muon and tau 
g — 2 measurements do not pose any bound on this scenario. 

Flavor physics constraints have also been studied in the literature. The charged scalars 
will mediate lepton flavor violating decays. In g —> ey, for instance, the additional branch¬ 
ing ratio is proportional to (m H ±v 2 )~ 4 - Because of that, it is always possible to evade this 
bound for large enough values of V 2 (or m^±). The limits we derive here are independent 
of a particular choice of V 2 as they concern directly the spectrum. 

4 Analysis of the Models 

For each model we generate ~ 10 7 points. For each of those points we calculate the 
corresponding scalar potential parameters and verify if they fulfill the constraints described 
in sec. 3. We only show on our plots the allowed points, which are about 10% of the 
generated ones. Unless stated otherwise, the points are color coded accordingly to the fit 
to EWPT data: blue, green, and red correspond to the 1, 2, and 3a allowed regions, while 
gray points are excluded at 3er or more. 

2HDM with a Z 2 symmetry. Let us first discuss the results for the 2HDM with a 
Z 2 symmetry. As discussed in sec. 2.1, the model has a very light neutral scalar. In fact, we 
verified that eq. (A.7) and the perturbative unitarity conditions (3.5) require mgj lOx-i^. 
Moreover, as the scalar potential parameters Aj and rri 1 - can be written in terms of the 
physical masses and the vevs, we perform a scan in the physical parameter space, imposing 

5 There would be a small contribution due to modifications of the h —> 77 coupling, but Higgs data 

already constraint it to the level that there is no observable modification to the muon g — 2 . 
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the following conditions 


0.01 eV <mu < 1 GeV, 

124.85 GeV <m h < 125.33 GeV, 

70 GeV <m H ± < 1 TeV, 

1 GeV <m,A < 1 TeV, 

— 7r/2 <a < it/2, 

0.01 eV <V 2 < 1 MeV. 

The Higgs mass range is taken from the ATLAS+CMS measurements combination in 
ref. [53]. Note that although a has to be small we did a scan over the whole physical 
range of this parameter, since we wanted to be as general as possible. By using a logarith¬ 
mic prior, for convenience, we found that indeed a is small. 

Using the power of perturbative unitarity constraints, we found that the CP-odd and 
charged scalars are restricted to be below ~ 600 — 700 GeV. This can be easily understood 
from eqs. (A. 8 ) and (A.9). Since M 2 oc m\ 2 = 0 and A 4,5 cannot be too large, the masses 
cannot go arbitrarily above the electroweak vev. 

Moreover, the presence of a very light scalar in the spectrum, below the GeV scale, 
yields a substantial negative contribution to the S parameter. The impact of the EWPT 
can be seen in fig. (1), where all points scanned were projected in the S x T plane and the 
allowed region by EWPT was drawn. Remarkably, only very few points (in red) were found 
which provide a viable model, within the 3a allowed region for the EWPT. From our scan, 
it can be concluded that: the T parameter strongly prefers tua ~ or a lighter H± 

with m H ± ~ 150 GeV together with a via > 300 GeV; while the S parameter, although it 
depends very mildly on the charged and pseudoscalar masses, exhibits a slight preference to 
this latter region. All in all, the values of S are always below ~ —0.25, revealing a tension 
with EWPT always above the 2.97a level . 6 As an example, we obtained the following 
scalar spectrum, which is allowed at 2.99a: 

mn = 0.18 eV, = 124.9 GeV, m H ± = 158 GeV, 
itla = 567 GeV, tana = —9.3 x 1CU 6 , tan/3 = 2.3 x 10 -6 . 

From this analysis, we can conclude that the 2HDM with a Z 2 symmetry is definitely very 
disfavored by data. It is not even clear that the region found which is in the 3a border 
of EWPT is really viable. A closer look into this region of the parameter space reveals 
that these points suffer from at least one of the following worrisome situations: (i) the e\ 
scattering amplitude, in eq. (3.5), is on the verge of violating unitarity, with at least about 
~ 98% of the bound saturated; (ii) the same for a + scattering amplitude, with at least 
~ 98% of the bound saturated; (iii) the stability condition is very fragile, with the third 
condition of eq. (3.1) satisfied with a relative difference of less than ~ 4 x 1CU 4 ; and (iv) 

6 To be precise about such strong statement, we also included in our analysis the accepted points of a 
second scan centered on the red region, where the charged scalar mass range was changed to 150-160 GeV 
and the pseudoscalar mass range was changed to 500-580 GeV, with 10 5 points. 
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Figure 2. Neutrinophilic 2HDM with Z 2 symmetry. Left: Predicted values for S and T (left) and 
isolines of S and T values as a function of rriA and m#±, for rriH mz- Right: tana x tan/3 
plane exclusions obtained using h invisible width. Orange points are excluded, while blue points 
are allowed. 

the same but for the second condition of eq. (3.1), satisfied with a relative difference of 
less than ~ 0.05. Therefore, given this delicate region of the parameter space, it would be 
important to include radiative corrections to see if the stability and unitarity of the model 
still holds at one loop. Notice that by using 47r as the perturbative unitarity limit this 
small region disappears. 

A possible way to evade these problems could be to have a larger so that the mass 
spectrum, specially m#, becomes more flexible. Nevertheless, unless u 2 > O(GeV), the 
problem does not disappear, strongly disfavoring this minimal model as an explanation for 
neutrino masses. 

One could now be tempted to include a right-handed neutrino contribution, dropping 
the lepton number conservation symmetry of the model. In fact, as u 2 is small, it may 
be possible to have a low-energy realization of the type I seesaw scenario which leads to 
observable sterile neutrino phenomenology, and hopefully could increase a bit the value of 
the S parameter to make the model viable. As the effect on S grows with the mass of 
the fermions in the loop, we make a distinction between two regimes: the right-handed 
neutrinos can be below or above the GeV scale. In the first, what happens is that the 
contribution to the S parameter is suppressed by the ratio between these small masses and 
the Z mass and can be neglected (for instance, the active neutrino contribution to S is 
virtually zero). In the second case, although the sterile neutrino masses might be large, 
the coupling to the Z is suppressed by the active-sterile mixing which generically goes as 
the ratio between the active to sterile neutrino masses, Therefore the impact of 

right-handed neutrinos is never large enough to substantially change the S parameter'. 

For completeness, we also show in the right panel of figure (2) the impact of the Higgs 
invisible width measurement in the tana x tan/3 plane. Given the preference for heavier 
S, we will consider the case where only h —> HH is present. From eq. (3.11), since m\ 2 = 0 

'This fact has also been checked numerically using the expressions in ref. [54]. 
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in this model, the QhHH coupling can be rewritten in the limit of small f3 and a as 


9hHH ~ 


m\ sin( 2 a) 
v sin(2/3) ’ 


(4.1) 


which can be sizable only if a > (5, explaining the behavior of the excluded region (orange) 
in fig. ( 2 ). Since the ratio a/(3 is already constrained by the theoretical limits (see eq. (2.9)), 
this constraint turns out to be less stringent than the others. As a last comment, the 
charged scalars could also have an impact on h —> 77. In the small 3 a allowed region, 
the modifications to the diphoton width are generically between ± 10 %, depending on the 
precise values of A 3 . This quartic coupling only affects rn h , so it is only weakly bounded 
by perturbative unitarity. 

Finally, since the smallness of mu causes the tension with EWPT, one may wonder 
what is the impact of loop corrections on the scalar spectrum of this model. Generically, 
in a 2HDM with Z 2 symmetry, the charged and CP-odd mass matrices are not modified by 
one-loop corrections. The CP-even matrix receives radiative corrections of the form [55] 


/ Ai^i A 3 45 nin 2 j 1 ( Am^vf Am\ 2 viv 2 \ 
l A 34 5 ^iT 2 A 2 U 2 ) 6471-2 l Am\ 2 v\v 2 Am 22 v 2 J ’ 


(4.2) 


where the second term comes from the one-loop effective potential. Since Am| are solely 
functions of masses and quartic couplings, the dependence of the CP-even mass matrix on 
W 72 is preserved at one-loop level, implying a small value for mu if v 2 is small. We have 
checked by explicit calculations that the corrections to mu are at the most a factor 100 , 
which is still insufficient to solve the problem with the S parameter. 

2HDM with a global U ( 1 ) symmetry. We now focus on the phenomenology of 
the softly broken 17(1) model. A non-zero m\ 2 term allows for heavier H, presenting a 
major change in the phenomenology with respect to the previous model. Without the 
requirement of a light scalar, we enlarge the scanned region accordingly. The absence of 
the A 5 quartic coupling makes the pseudoscalar degenerate in mass with H (to first order 
in v 2 ). Therefore we perform an initial scan of the spectrum parameter space, this time in 
the region 


10 GeV <mu < 1 TeV, 
124.85 GeV <m h < 125.33 GeV, 
70 GeV <m H ± < 1 TeV, 
m A = m Hl 
—7 t/2 <a < n/2, 

0.01 eV <v 2 < 1 MeV, 


as well as a second scan with m#± and m A heavier then 1 TeV and almost degenerate. 
We follow the same procedure as before, showing only the points allowed by perturbative 
unitarity and stability constraints. The results are presented in fig. (3). 

In contrast to the previous case, due to the possibility of obtaining a heavier H in the 
mass spectrum, there is a region of the parameter space of this model which passes the 
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Figure 3. Neutrinopliilic 2HDM with softly broken global 17(1) symmetry. The blue, green and 
red points are allowed by EWPT at la, 2cr, and 3cr, respectively, while the gray points are ruled 
out at 3cr. Top left: parameter space in the plane mu x m#± which satisfy perturbativity, unitarity 
and stability constraints. Top right: projection of these points in the S x T plane. Bottom left: 
h —?> 77 signal strength as a function of mu±- Bottom right: region in the mu x v 2 plane that is 
excluded by the Z invisible width (orange points). 

electroweak precision tests and theoretical constraints. The behavior of the T parameter 
is similar to the previous scenario: either the mass splitting between A and H± is at most 
~ 80 GeV, or the charged scalar is around 100 GeV while mu = rriA > 150 GeV, with 
negative values of T for larger mjj±. This explains the strong correlation on the allowed 
region in the upper left panel of fig. (3). We also present the projection of these points in 
the S x T plane in the upper right panel of fig. (3). For the a and (3 parameters we find 
that the allowed region is tan (3 < 1CP 6 and a < 5/3. 

As discussed in the previous sections, this model can also accommodate a pair of 
neutral scalars ( S = H,A ) satisfying ms < mh/2 if m\ 2 is small enough. In this case, 
the constraints coming from the Higgs invisible decays are similar to those described for 
the model with a Z 2 symmetry and turn out to be relatively weak. On other hand, the Z 
invisible width can provide valuable constraints when the channel Z —> Svu is open. To 
perform this analysis we scan over the oscillation parameters, imposing the perturbativity 
condition < ms/2. We show on the bottom right panel of figure (3) the excluded region 
(orange points) under these assumptions in the mu X V 2 plane. The region ms < mz/2 
is completely excluded, because in this case we integrate over the poles of the off-shell 
scalars in Z -A H{A* -A isP) and Z —> A(H* -A z/P), enhancing the decay rate by orders of 
magnitude. 

For a heavy enough , as can be seen in the lower left panel of fig. (3), the h —> 77 
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signal strength is diminished by about ~ 5%. We can understand this non decoupling 
feature by noticing that the hH + H~ coupling is —iX^v, which in turn has a correlation 
with rn,fj±, specially in the larger mass region. This can be understood by noticing that, 
in eq. (A.3), for large m H ±, we have 


A3 


sin 2a \ m 2 H ± 
sin 2/3 J v 2 


(4.3) 


Typically, a < 5/3, which corresponds to a strong correlation between A 3 and m H ±, and 
this is the denser region around // 77 = 0.95. However this is not always the case, and the 
correlation is lost when the ratio of sines is closer to 1 , now corresponding to the sparser 
points with a much weaker correlation. Nevertheless, we see that for a heavy enough 
charged scalar, the contribution to the Higgs diphoton width is always negative. 

One could ask if it is also possible to have a Majorana mass term, since the 17(1) 
symmetry is softly broken. First, as pointed out in refs. [54, 56], the impact of heavy 
right-handed neutrinos via loop effects on electroweak precision observables is very small. 
Therefore, there is no significant interplay between this and the scalar sector of the model, 
and thus the phenomenology studied here would be essentially unchanged. On the other 
hand, if we consider an UV completion that simultaneously breaks the symmetry and 
originates a Majorana mass term, we find that such scenario is non-minimal, i.e., at least 
two new fields have to be included. 


Comments on non-minimal models 

Due to the large number of possible variants, performing exhaustive analises of non- 
minimial models is unpractical and well beyond the purpose of this paper. Nonetheless, 
we may glimpse the phenomenology of some representative cases. 

A neutrinophilic 2HDM with a softly broken Z 2 symmetry would surely be allowed 
by data, in contrast to the spontaneously broken Z 2 scenario. Such model would be more 
general than the two models considered here, as it would span a non-zero value of both 
mi 2 and A 5 . As can be seen from eqs. (A.7) and (A. 8 ), the simultaneous presence of these 
terms in the scalar potential lifts the degeneracy between the neutral scalars A and H. In 
fact, we have checked that the allowed region in the plane tiia x mjy± is very similar to 
the one exhibited in fig. (3) (top left panel), except for the fact that the T parameter now 
implies a correlation only between m^± and rriA- If the neutral scalar H decays dominantly 
to neutrinos, as it is likely to happen, it would be very difficult to probe it by resonant 
production at colliders. 

Another way of evading our limits would be to enlarge the particle spectrum of the 
spontaneously broken Z 2 model (or generically any Zjv), for instance, by adding a scalar 
singlet S, doublet T 3 or triplet A, all charged under the new symmetry. In the singlet 
case, a triple or quartic term S' < hj < I >2 or S' 2 < h{ < I > 2 could be present in the potential for 
a judicious choice of charges. After the singlet acquires a vev, this term would play a 
role similar to the 77112 soft breaking term allowing for larger values of mn ■ However, 
the quartic always present, would induce a Higgs-singlet mixing. This would 

diminish all Higgs couplings to fermions and gauge bosons by a factor sin# where 9 is the 
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corresponding mixing angle. The mixing is constrained by Higgs production cross section 
measurements to be sin 2 6 < 0.2 [45]. In the case of adding a third doublet, the triple 
coupling is impossible, but a quartic one could be present. Last, in the case of a scalar 
triplet, although a triple coupling would be possible, a large triplet vev could irrevocably 
disturb electroweak precision tests, especially the T parameter. 

5 Conclusion 

We performed an analysis of the minimal neutrinophilic two-Higgs-doublet models which 
can accommodate neutrino masses by means of the tiny vev of the additional Higgs dou¬ 
blet. The models studied here differ among themselves by the symmetry that forbids the 
couplings between neutrinos and the scalar which gets the electroweak scale vev. The cases 
studied here span a discrete Z 2 and a softly broken global 17(1) symmetry. 

The bounds considered come both from theory and experiment. The unitarity per¬ 
turbative requirement at tree-level strongly constrains the scalar mass spectrum of these 
models, either by the presence of a very light neutral scalar {mu ~ V 2 ), in the Z 2 model, 
or with a degeneracy between the scalar and pseudoscalar particle masses (mu = ttia), in 
the global 17(1) scenario. 

If there is no additional particle content, the Z 2 symmetry model was found to be in se¬ 
vere tension with the electroweak precision tests, due to the very light neutral scalar, which 
generates a large negative contribution to the S parameter. The inclusion of a Majorana 
mass term for the right-handed neutrinos, providing a low scale realization of the seesaw 
type I mechanism, does not save the model, as the right-handed neutrino contribution to 
the S parameter is always negligible. Therefore, we conclude that the neutrinophilic 2HDM 
with a spontaneously broken Z 2 symmetry is strongly disfavored by data. 

The analysis of the model with an explicit broken global 17(1) symmetry reveals a 
region of the parameter space which is allowed by all bounds considered. Due to the 
set of constraints and the symmetries of the model itself, the spectrum is quite limited. 
The U( 1) symmetry predicts that the neutrinophilic scalar is degenerate in mass with the 
pseudoscalar, mu = mA- Besides, the electroweak precision tests play a very important 
role, specially the T parameter which is sensitive to the absolute mass splitting of the 
pseudoscalar and the charged scalars, limiting it to be at most ~ 80 GeV. Therefore, 
an important consequence of the theoretical and experimental constraints is that, if the 
new scalars are above ~ 400 GeV, all these particles should have very similar masses. 
Moreover, the Z invisible width excludes the region mu = mA < mz/ 2. Besides, the 
h —> 77 branching fraction might be modified by about ±30% for m#± < 200 GeV, while 
for heavier above 500 GeV, this ratio can be atmost 1 or lower by 5%. Finally, we 
stress that this model can be well within the reach of LHC 13 TeV, by probing the h —> 77 
branching fraction of by direct pair production of the charged scalars, if they are below 
0(300 GeV). 
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A Tree-level Relations for the Quartic Couplings 


The quartic couplings can be expressed in terms of the physical masses, vevs and mixing 
angles as: [24] 


Ai 

A 2 

A3 

A 4 

A 5 


\ ( — tan 2 (3M Z + 


1 


sin 2 a 


cos 2 (3 


m H + 


cos 2 a 


9 /a it r 1 ) COS Ot O 

— cot 2 (3M Z H— 9 n m H + 


cos 2 (3 
sin 2 a 


ml 


—M 2 T 2m~(j± T 


sin 2 (3 
sin( 2 a) 


m l 


sin (2 p) 


sin 2 f3 
(m 2 h - m \) 


^2 (M 2 + m\- 2m 2 H± ) 

^2 (M 2 - mi) , 


where M 2 


m 12 

sin f3 cos (3 


Inversely, we have 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 


ml = M 2 sin 2 (a — (3) 


+ ( Ai cos 2 a cos 2 /3 + X 2 sin 2 a sin 2 (3 + sin 2a sin 2/3 ] v 2 , (A. 6 ) 

mjj = M z cos 2 (a — (3) 

+ (Ai sin 2 a cos 2 (3 + X 2 cos 2 a sin 2 (3 — —^ sin 2 a sin 2(3) v 2 , 


m\ = M 2 — A 5 U 2 , 

2 1, f 2 A45 2 

m H± = M - —v , 


where A 45 = A 4 + A 5 and A 345 = A 3 + A 4 + A 5 . 


(A.7) 

(A. 8 ) 

(A.9) 
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